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Abstract 

We present an emergent universe scenario making use of a new solution of the 
Starobinsky model. The solution belongs to a one parameter family of solutions, where 
the parameter is determined by the number and the species (spin- values) of primordial 
fields. The general features of the model have also been studied. 
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1 INTRODUCTION 



Recently Ellis and Maartens [1] reconsidered the possibility of a cosmological model in 
which there is no big-bang singularity, no beginning of time and the universe effectively 
avoids a quantum regime for space-time by staying large at all times. The universe starts 
out in the infinite past as an almost static universe and expands slowly, eventually evolving 
into a hot big-bang era. An interesting example of this scenario is given by Ellis, Murugan 
and Tsagas [2], for a closed universe model with a minimally coupled scalar field </>, which 
has a special form of interaction potential V{4>). It was pointed out that this potential is 
similar to what one obtains from a R+aR 2 theory after a suitable conformal transformation 
and identifying <p = — a/3 ln(l + 2aR) with a negative a. Although the probability of these 
solutions is not high, the emergent universe scenario nevertheless merits attention as it 
solves many conceptual and technical problems of the big bang model. In this paper, we 
point out that the Starobinsky model, the original as well as the modified version, permit 
solutions describing an emergent universe. The solution may be used to model varieties 
of cosmological scenarios consistent with the observational results, available at present 
and expected in the near future. Thus, it may be possible to build models which avoid 
the quantum regime for space-time but share the good features of the standard big bang 
model. 

In the models considered in Ref . [1] and [2], a closed universe was considered. However, 
recent results from BOOMERANG and WMAP indicate that the universe is most likely to 
be spatially flat. If the universe has always been large enough, the field equations become 
simpler. In fact, in the Starobinsky model, the field equations can be written as a second 
order differential equation for the Hubble parameter H, vide equation (7). 

2 STAROBINSKY MODEL : 

In the Starobinsky model, one considers the semi-classical Einstein equation, 

R»v - \g^R = -8itG< Tp, > (l) 

where < T^ u > is the vacuum expectation value of the energy momentum tensor. In the 

case of free, massless, conformally invariant fields the vacuum expectation value < T^ v > 
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can be written, with the Robertson Walker metric as, 



(2) 



where K\ and K3 are numbers and 





(3) 
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Equation (2), written for < >, gives the well-known trace-anomaly, indicating that the 
conformal invariance is broken by the regularization process. Note that the two tensors 
^H^ u and ^H^ u have different features. The tensor is identically conserved and 

can be obtained by varying a local action A ~ / \f^g R 2 d 4 x. However, one of the counter 
terms to be added to the lagrangian to regularise < > has the form Hyf^g R 2 , where 
H is a logarithmically divergent constant, thus permitting the possibility of an addition 
of any finite number to fx. Hence, K\ is not determined and can be given any arbitrary 
value including a negative one. The tensor is conserved only in a conformally flat 

universe and it cannot be obtained by variation of any local action. The constant K% is 
fully determined, once the fields are specified, e.g., 



where Nj gives the number of quantum fields of spin I. The freedom of choosing K\ 
arbitrarily has led to the modified Starobinsky [3,4] model, where one adds a counter term 
~ \f~~9 R 2 , and chooses K\ » K3. The theory then becomes essentially a typical R 2 
-theory. However, the earlier work on Starobinsky models were done in the context of a 
big bang model. The emergent universe model needs a different scenario and one makes 
different choices of K\ and K3, as will be shown below. 

Let us choose 8itG = 1 and write the evolution equation (1) of the flat FRW universe 

as 




(5) 



which can be written in terms of the Hubble parameter 

H 2(JL- H A = -^1(2HH + 6H 2 H-H 2 ). (7) 

Equation (7) will determine the evolution of the universe when the initial conditions are 
provided. 

3 EMERGENT UNIVERSE : 

We now look for a solution which describes a universe which exists eternally and always 
remains large so that a classical description is possible at all times. The matter, as indicated 
earlier, will be described by quantum field theories. We first note the following two exact 
solutions : 

• H = is a solution but it describes an eternally static universe. It can be checked 
that the solution is stable against small linear perturbations, if K\ is positive. However, 
with K\ negative, the state is unstable. 

• H = gives a de Sitter type solution, but the solution is not stable. 
To get an emergent universe, we look for a solution of the suggestive form 

a(t) = a ((3 + e at Y (8) 



where the constants a and uj will be determined from equation (7). The constants a Q and 
(5 may be determined from initial conditions. Since we have 

_ uae at ■ _ uja 2 e at » _ u;a 3 e at (l - e at ) 

~ f3 + e at ' ~ (/3 + e«*) 2 ' ~ W + e at f W 

the function in (8) will be a solution, if u = |, a = §^/5j and K\ = —-^K^. Thus K\ 
will be chosen negative in this model. 

The fact that the emergent universe scenario is indeed permitted by equation (7) comes 
as a surprise. Since K\ and K% are related, we have a one parameter family of solutions 
specified by the value of K%. The general features of the solutions are 

1) The scale factor a(t) has a non-zero value a(t — ► — oo) = a /? 2 / 3 as t — > — oo, which 

may be chosen much larger than the Planck length. The Hubble parameter H and its 
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Figure 1: Variation of ft. = HX10 2 and the derivatives, H and H are represented by thin, 
broken and thick lines respectively ( (3 = 10000 and = 1.8) 



time derivatives H and H all vanish in the limit t — ► — oo. Thus the solution decribes an 
emergent universe with natural initial conditions. 

2) If (3 » 1, the universe remains almost static during the period — oo < t < 
2 /X3 ln/3 = t Q . Thus during the entire infinite past of t = 0, the universe expands 
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only by a factor 



a(t = 0) ( 1\ 2/3 

—-i '— = 1 + - ~ 1 (10) 

a(*->-oo) \ PJ 



for a large /3. 



For i > |V^3 In/?, the universe expands rapidly, eventually reaching an asymptotically 
de Sitter stage. Fig 1, gives the time variation of the Hubble parameter and its derivatives, 
for K3 = 1.8 (which occurs if the particle number and species correspond to the minimal 
SU(5) model) and /? is given a value 10 4 . 

3) It may be useful to check if the solution is stable under small linear perturbations. 
We write 

H = H 8 (1 + S) (11) 
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where S is a small perturbation and H s = ^ +eat ^ is the solution obtained above. Substi- 
tuting this in equation (7) , we obtain the equation for 5 : 

S + A(t)6 + B(t)5 = (12) 

where 

Note that § < A(t) < 2a and B(t) > for -oo < t < oo. It is difficult to solve the 
equation (12). However, if we assume that 8 has a solution 5 ~ e m *, we must have, for t 
negative, 

m 2 + %n ~ (14) 

and m cannot be positive, since a and /3 are both positive. Thus the solution seems to be 
stable under small perturbations at least in the negative t regime. 

4) The evolution can be described alternatively in terms of a scalar field by substituting 
H = cj? in equation (7). This gives 

+ 3^ + ^ = (15) 

d(p 

where V = —m?4> 2 + \<j) 6 with m 2 = A8 \~ Kl \ and A = l4 j^ 3 Kl | ■ Note that V{<p) has the shape 
of a double well potential with a maximum at (f> = and two minima for <p m = ±(^) 1 / 4 , 
giving V m i n = — 72 \xI\^/k 3 ' ^ so no * e ^ = ^ gi yes an unstable static universe while the 
evolution of the emergent universe is given by the path from V{4>) = to V m in- However, 
H here is time dependent and is determined self-consistently by the instantaneous value of 
<p. The equation (15) does not make the problem simpler, although it looks very familiar 
and attractive. 

5) The prescribed evolution of the universe gets modified when particles are produced 
due to the expansion of the universe. We have considered above only massless, conformally 
coupled fields in conformally flat spacetimes, i.e., the conformally trivial situation where 



no particles are produced. The perturbative calculations of particle production make use 
of small deviations from this conformal triviality. One considers either (i) a small mass 
and/or (ii) deviations from a conformally flat space time, say by assuming a Bianchi I 
spacetime deviating a little from a FRW spacetime or (iii) assuming a small non conformal 
coupling with |£ — || = e. For simplicity, we shall consider here the case of non conformal 
coupling. Let us consider a real massive scalar field satisfying the equation 

<j> + 3H<j> + (m 2 + £R)<f> = (16) 

with |£ — g| = e << 1. As the scale factor changes, the changing gravitational field feeds 
energy into the prturbed scalar field modes. As long as the mode frequency of the field 
is greater than the Hubble expansion rate, a co-moving detector will not respond. Modes 
with lower frequencies will, however, be excited. Thus the presence of mass makes the 
particle production process less efficient. To begin with we consider the case m = and 
e ^ 0. Zeldovich and Starobinsky, Starobinsky [3], Birrell and Davies [3] and Vilenkin [4] 
have studied this case. The rate of particle production per unit volume per unit conformal 
time n, (dn = ^fy) , is given by 

P = (17) 
dn 2 y J 

where R is the scalar curvature. This gives 



2a n dn 



Y = 

e 2 dt 



36u 2 a 4 e 2at {l + 2Lue at ) 2 
(J3 + e «*) 14 /3 



(18) 



In fig. 2, we have plotted Y against t for /3 = 10, 000 and K3 = 1.8 and 2.5. We note that 

the rate of particle production peaks around t ~ t Q , as expected. Increase in the number 

of species of particles i.e., a higher value of K% will indicate a delay in particle production 

and also a wider peak. The dependence of the rate of particle production on the parameter 

(3 is shown in fig. 3. The absolute value of the peaks depends on the value of the scale 

factor a(t) and the constant a Q may be determined from the rate of particle production. 

Although we have considered zero mass particle we expect similar results even when m 

is nonzero but small ( see Vilenkin [ 4 ]). During the subsequent evolution the particles 

produced thermalize and the universe enters eventually into a radiation dominated stage 
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Figure 2: Dependence of rate of particle production on time for f3 = 10000 and K% = 2.5 
(broken line) and K3 = 1.8 ( thin line ) 

of the standard hot bigbang model. The details of the process of particle production and 
the evolution of the universe will be considered elsewhere. 

4 DISCUSSION 

We have presented here a one parameter family of solutions of the Starobinsky model which 
describes an emergent universe. The earlier work on Starobinsky model was done in the 
context of big bang singularity and the particle production was achieved by its oscillatory 
solution. The present solution on the other hand describes the universe which is almost 
dormant during the infinite past period — 00 < t < § hi/3, after which it undergoes a 
rapid expansion, see fig. 1. The parameters a Q and (3 are related with initial conditions. 
They also introduce a new mass scale in the process of particle production. We have 
not considered here any specific theory of particle interactions which will determine the 
subsequent evolution of the universe as well as the details of large scale structure formation. 
Whether the present scenario can successfully explain the present observational data needs 
further study. However, the fact that one encounters solutions describing an emergent 
universe in different contexts may be a good reason for taking such solutions seriously, 
although the probabilities for such solutions may not be very high. 

8 



Figure 3: Dependence of rate of particle production on time for K% = 1.8 and f3 = 10000 
(thin line) and (5 = 20, 000 ( broken line) 
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